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I.  INTRODUCTION 


Because  IMPATT  diodes  operate  at  Junction  temperatures  well  in  excess 
of  the  operating  temperatures  customary  for  other  semiconductor  devices, 
they  frequently  fail  prematurely  as  a  result  of  unexpected  thermal 
problems.  Therefore,  it  is  important  to  model  IMPATT  diodes  thermally  in 
order  to  be  able  to  predict  the  life  expectancy  of  a  device.  For  this  an 
accurate  thermal  model  is  needed,  because  the  normal  operating  temperatures 
of  IMPATT  diodes  frequently  approach  critical  levels.  Although  no  model 
can  replicate  operating  conditions  exactly,  the  best  models  approximate  the 
true  situation  accurately  enough  so  that  the  conclusions  drawn  by  the 
modeler  are  valid. 

As  is  well  known  from  elementary  physics,  when  different  parts  of  a 
body  are  at  different  temperatures,  heat  flows  from  the  hotter  parts  toward 
the  colder  ones.  In  the  case  of  IMPATT  diodes,  heat  is  generated  at  the 
p/n  Junction,  so  heat  flows  from  the  junction  toward  the  heat  sink.  This 
flow  of  heat  can  occur  in  several  different  ways,  but  in  the  case  of  IMPATT 
diodes  the  most  important  method  of  heat  transfer  is  conduction.  By  defi¬ 
nition,  heat  conduction  is  the  transfer  method  by  which  heat  passes  through 
the  body  itself.  The  model  developed  for  double-drift  IMPATT  diodes  limits 
the  heat  flow  to  heat  conduction  and  neglects  all  other  heat  transfer 
methods,  such  as  radiation  and  so  on.  By  limiting  heat  flow  to  conduction, 
one  can  develop  a  thermal  model  based  on  the  thermal  resistance.  By  defi¬ 
nition,  the  thermal  resistance  6  of  any  structure  can  be  obtained  from  the 
relationship 


6  =  AT/P  [K/W]  (1) 

where  AT  is  the  rise  of  temperature  due  to  the  dissipated  power  P,  which  is 
converted  into  heat.  Note  that  e  is  the  inverse  of  the  thermal  con¬ 
ductance. 
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In  this  report  we  develop  a  thermal  model,  based  on  the  thermal 
resistance  of  the  IMPATT  structure,  that  is  applicable  to  double-drift 
IMPATT  diodes.  However,  prior  to  the  development  of  the  thermal  model,  it 
is  advantageous  to  review  briefly  some  important  features  of  IMPATT  diodes 
and  the  diode  structure  itself. 
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II.  IMPATT  DIODE  STRUCTURE 


In  1958  W.  T.  Read1  predicted  that  a  localized  high-fiela  avalanche 
region  followed  by  a  moderately  high-field  drift  region  would  make  an  effi¬ 
cient  source  of  microwave  energy.  It  took  nearly  ten  years  after  Read's 
paper  for  silicon  (Si)  technology  to  mature  to  the  point  that  a  Read  diode 
structure  could  be  fabricated.  The  results  showed  that  Read's  original 
predictions  were  correct,  but  they  also  showed  that  the  elaborate  outline 
proposed  by  Read  was  unnecessary.  Since  then,  a  variety  of  devices  based 
on  Read's  principle  for  generating  microwave  power  have  been  developed, 
most  recently  the  double-drift  IMPATT  diode. 

A  double-drift  IMPATT  diode  consists  of  a  heavily  doped  p4  region 
followed  by  a  moderately  doped  p  region,  a  moderately  doped  n  region,  and  a 
heavily  doped  n4  region.  Avalanche  breakdown  occurs  at  the  p/n  Junction, 
while  the  n4  and  p4  regions  allow  ohmic  contacts  to  be  made  to  the  IMPATT 
diode.  Figure  1  shows  the  doping  profile  of  a  p4/p/n/n4  double-drift 
IMPATT  diode,  along  with  the  electric  field  at  avalanche  breakdown.2 

During  the  operation  of  an  IMPATT  diode  only  a  relatively  small 
portion  of  the  applied  dc  power  is  converted  into  microwave  energy;  the 
rest  is  converted  into  heat.  The  ratio  of  the  microwave  energy  to  the 
total  applied  dc  power  is  called  the  conversion  efficiency  of  the  IMPATT 
diode.  To  obtain  acceptable  conversion  efficiencies,  high  current  densi¬ 
ties  are  required  under  avalanche  operating  conditions;  i.e.,  efficient 
operation  requires  large  power  dissipation.  Under  continuous  operation 
this  power  dissipation  can  be  achieved  only  via  structures  that  are  capable 
of  handling  the  power. 

In  practical  applications  there  is  a  specified  maximum  operating 
temperature,  dictated  by  reliability  considerations,  that  corresponds  to  a 
AT  ;  at  the  same  time,  there  is  a  desired  conversion  efficiency  that 
dictates  a  minimum  current  density  Jmin.  Because  of  the  AT  x  and  the  Jmin 
requirements,  continuous  operation  can  be  achieved  only  by  small  diode 
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chips.  Therefore,  IMPATT  diodes  are  fabricated  by  means  of  a  mesa  outline 
having  a  small  chip  radius  of  a  <30  ym.  An  equation  defining  the  maximum 
chip  radius,  a,^,  relative  to  ATmax  and  is  given  later. 

To  reduce  the  thermal  resistance  of  the  assembly,  the  IMPATT  diode  is 
bonded  with  its  Junction  side  down  to  a  good  thermal  conductor;  this 
reduces  the  distance  the  heat  flux  must  travel  from  the  heat-generating, 
avalanching  p/n  Junction  to  the  heat  sink.  In  the  case  of  double-drift 
IMPATT  diodes,  the  heat  sink  consists  of  a  small  cylindrical  piece  of 
diamond  (Type  II)  that  is  hot  pressed  into  a  copper  (Cu)  base.  Diamond  is 
Used  for  its  excellent  heat  conductivity,  the  highest  known  to  date  in  the 
temperature  range  where  IMPATT  diodes  operate.  Thinning  the  n+  substrate 
to  <25  um  not  only  reduces  the  series  ohmic  resistance  but  also  improves 
the  conversion  efficiency  of  the  IMPATT  diode. 
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II.  DIODE  THERMAL  MODEL 


The  thermal  resistance  of  a  single-drift  IMPATT  structure  has  been 
analyzed  before.  Swan  et  al.^  have  shown  that  the  total  thermal  resis¬ 
tance  eT  of  a  single-drift  IMPATT  structure  can  be  expressed  as  the  sum  of 
the  thermal  resistances  of  the  individual  components,  provided  there  is  an 
ideal  thermal  contact  between  all  interfaces.  Their  analysis  is  general 
and  therefore  applicable  to  double-drift  configurations. 


Figure  2  shows  a  simplified  outline  of  a  double-drift  Si  IMPATT 
structure  mounted  on  a  diamond  heat  sink;  this  outline  is  used  to  develop 
the  thermal  model.  Note  that  the  dimensional  relationships  are  distorted. 
Figure  2  shows  that  there  are  two  heat-conducting  paths,  (1)  downward 
toward  the  heat  sink  and  (2)  upward  through  the  top  contact.  Therefore, 
the  total  thermal  resistance  of  the  IMPATT  structure  is  the  sum  of  these 
two  and  can  be  expressed  as 


=  e  .  ♦  e 

down  up 


edo «, 


(2) 


The  thermal  resistance  of  the  upward  path  has  been  neglected,  since  it 
provides  less  than  51  of  the  total  heat  conduction.  Thus  the  total  thermal 
resistance  of  the  IMPATT  diode,  6^,  can  be  expressed  by  the  thermal 
resistance  toward  the  heat  sink  as  the  sum  of  four  components: 


eT  *  eSi  *  emetal  *  edia  *  eCu  lK/“! 

w h*r  >  Bgj  is  the  thermal  resistance  of  the  mesa  diode,  C^tai  is  the 
thermal  resistance  of  the  metal  layer,  edia  is  the  thermal  resistance  of 
the  diamond  heat  sink,  and  eCu  is  the  thermal  resistance  of  the  copper 
base.  Consequently,  we  simplify  the  development  of  the  thermal  model  by 
defining  only  the  components  of  Eq.  (3). 
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Fig.  2.  Diagram  of  am  IMPATT  Diode  Chip  Mounted  on  a  Diamond  Heat  Sink 


Gibbons  and  Misawa*1  have  analyzed  the  temperature  and  current  distri¬ 
bution  of  an  avalanching  p/n  Junction  on  a  homogeneous  semi-infinite  heat 
sink  having  thermal  conductivity  k.  They  used  a  model  in  which  the  diode 
was  represented  by  a  heat  flux  over  a  small  disc  of  radius  a;  the  heat  flux 
was  incident  upon  the  center  of  a  large  cylindrical  Cu  heat  sink.  They 
further  assumed  adiabatic  conditions  on  the  remainder  of  the  top  surface. 
The  heat  flux  was  the  product  of  the  junction  voltage  V,  which  was  assumed 
to  be  uniform,  and  the  current  density  J  (A/cm^),  which  is  a  function  of 
temperature.  They  concluded  that  as  the  input  power  is  increased,  the 
radial  distribution  of  the  current  density  becomes  less  uniform  while  the 
radial  temperature  distribution  becomes  more  uniform.  A  current  density  of 
J^in  results  in  an  approximately  uniform  heat  flux  over  the  disc.  Thus  the 
thermal  resistance  of  the  Si  mesa  diode,  6g^,  can  be  represented  by  a 
one-dimensional  heat-flow  model  as 


e 


si  *  V<kSi'a2> 


(4) 


where  lav  defines  the  distance  between  the  p/n  junction  and  the  heat  sink, 
kSi  is  the  thermal  conductivity  of  Si  (a  subject  we  will  discuss  later), 
and  a  is  the  radius  of  the  mesa  diode.  Their  analysis  is  general  and  is 
applicable  to  double-drift  designs  that  use  diamond  heat  sinks. 

Further,  it  is  assumed  that  spreading  of  the  heat  flux  in  the  thin 
metal  layers  of  the  diamond  heat  sink  is  negligible.  As  a  result,  when  one 
calculates  the  thermal  resistance  of  the  thin  metal  layers  of  the  diode 
chip  and  the  diamond  heat  sink,  the  two  resistances  can  be  combined  and 
represented  by  the  following  one-dimensional  heat-flow  model: 

W  *  ‘"“M,  Yki  (5> 

where  lj  and  k^  represent  the  thickness  and  thermal  conductivity,  respec¬ 
tively,  of  the  various  metal  layers;  n,  the  number  of  different  metal 
layers,  is  usually  three,  with  chromium  (Cr),  platinum  (Pt),  and  gold  (Au) 
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constituting  the  three  components.  These  metal  layers  are  applied  to  both 
the  Si  chip  (top  and  bottom)  and  the  diamond  heat  sink;  their  thicknesses 
are  given  in  Table  1 . 


Table  1.  The  Thickness  of  Metal ization  Layers  Applied 
to  Silicon  and  Diamond  Heat  Sink 


Metal 

Thermal 
Conductivity, 
ki(W/cm  K) 

Thickness , 
lj  (urn) 

Si  Chip 

Cr 

0.37 

0.06 

Pt 

0.73 

0.20 

Au 

3.15 

0.70 

Diamond  Heat  Sink 

Cr 

0.37 

0.06 

Pt 

0.73 

0.035 

Au 

3.15 

1.00 

In  a  typical  assembly  the  IMPATT  die  is  thermocompression  bonded  on 
top  of  the  diamond  heat  sink,  which  in  turn  has  been  hot  pressed  into  the 
Cu  base.  The  radius  of  the  diamond  heat  sink  R  and  its  height  D  are  at 
least  ten  times  larger  than  the  radius  of  the  diode  chip  a.  When  the 
diamond  heat  sink  is  hot  pressed  into  the  Cu  base,  the  operation  leaves  a 
gap  between  the  copper  and  the  side  of  the  diamond  heat  sink,  as  shown  in 
Fig.  2.  The  resulting  gap  is  small,  but  it  causes  heat  transfer  to  occur 
only  over  the  base  of  the  diamond,  not  over  its  circumferential  face. 

To  obtain  the  heat  sink  assembly  must  be  considered.  The  diode 
is  represented  by  a  uniform  heat  flux  that  enters  the  top  of  the  diamond 
heat  sink  over  a  small  centrally  located  area.  The  radius  of  the  heat  flux 
is  equal  to  the  mesa  radius  a.  The  heat  flux  can  exit  the  cylinder  only  at 
its  other  end,  since  the  remaining  surfaces  are  assumed  to  be  insulated. 
This  problem  has  been  analyzed  before. ^  Because  the  radius  R  and  the 
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height  D  of  the  diamond  heat  sink  are  much  larger  than  the  radius  of  the 
diode  chip  (i.e.,  R  >>  a  and  D  >>  a),  the  use  of  thermal  spreading 
resistance  is  Justified  when  expressing  e,jia.  Thus 

6dia  =  n/kdia,)  [(1/a)  *  (1/R)]  (6) 

where  k^ia  is  the  thermal  conductivity  of  diamond  (a  subject  to  be 
discussed  later),  and  a  and  R  are  the  radius  of  the  diode  chip  and  the 
diamond  heat  sink,  respectively. 

The  final  term  of  Eq.  (3)  is  the  thermal  resistance  of  the  Cu  base. 
Assuming  that  the  Cu  base  is  large  compared  to  the  radius  of  the  diamond 
cylinder,  we  can  express  6qu  by  using  the  thermal  spreading  resistance  of 
an  infinite  half-space.  Thus 


6Cu  *  1/<'*BkCu) 


(7) 


where  kCu  is  the  thermal  conductivity  of  Cu  and  R  is  the  radius  of  the 
cylindrical  diamond  heat  sink.  The  model  that  leads  to  Eq.  (7)  assumes 
that  the  temperature  at  the  Cu  and  diamond  interface  is  uniform.  This 
assumption  is  Justified  by  the  solution  given  by  Kennedy^  and  the  excellent 
thermal  conductivity  of  diamond  (Type  II). 
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IV.  THERMAL  CONDUCTIVITY  AND  THERMAL  RESISTANCE  OF  SILICON  AND  DIAMOND 


Equations  (4)  and  (6)  assume  that  the  thermal  conductivity  of  Si  and 
diamond  is  constant.  However,  the  thermal  conductivity  of  many  substances, 
including  Si  and  diamond,  is  not  constant  but  varies  as  a  function  of 
temperature.  In  the  temperature  range  of  interest  (i.e.,  the  normal 
operating  range  of  IMPATT  diodes,  from  room  temperature  to  a  few  hundred 
degrees  Kelvin  above  room  temperature),  the  thermal  conductivity  of  Si  and 
diamond  is  an  inverse  function  of  temperature: 

kSi  =  f(1/TSi)  (8) 


and 

kdia  =  f(1/Tdia>  <*> 

Figure  3  shows  the  thermal  conductivity,  based  on  actual  measured 
data,  of  Si,  Cu,  and  diamond  (Type  II). ^  Note  that  the  thermal 
conductivity  of  Cu  (as  well  as  all  metals)  is  practically  constant  in  the 
temperature  range  of  interest. 

In  calculating  the  thermal  resistance  of  the  Si  mesa  diode  and  the 
diamond  (Type  II)  heat  sink,  we  define  the  thermal  conductivity  of  Si  and 
diamond  by  the  following  empirical  formulas: 

kSi  =  292/ (TSi  -  Hi!)  =  292/TSe  (10) 


and 


kdia  =  4000/(Tdia  -  100)  =  H000/Tde  (11) 

where  TSi  and  Tdia  represent  the  temperatures  of  Si  and  diamond  in  degrees 
K,  respectively,  while  T^  and  Tde  express  effective  temperatures. 
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Fig.  3.  Thermal  Conductivity 


.  Equations  (10)  and  (11)  were  obtained  by  fitting  the  actual  curves  shown  in 

Fig.  3  at  two  points  each.  Comparative  plots  of  the  empirical  expressions 
of  the  thermal  conductivity  of  Si  and  diamond  (Type  II)  versus  the  actual 
measured  values  are  shown  in  Fig.  *4. 

Equations  (10)  and  (11)  relate  the  thermal  conductivity  of  Si  and 
diamond  to  the  actual  temperatures  inside  the  Si  chip  and  the  diamond  heat 
sink,  respectively.  Unfortunately,  during  operation  the  temperature  inside 
the  Si  chip  and  the  diamond  heat  sink  is  not  constant;  to  obtain  the  true 
value  of  eSi  and  t>dia,  the  temperature  distribution  of  the  Si  chip  and  the 
diamond  heat  sink  must  be  known.  As  this  would  be  difficult,  the  use  of  an 
approximation  simplifies  matters  greatly.  In  this  approximation  it  is 
assumed  that  fig^^  and  edia  can  be  represented  by  a  thermal  resistance  that 
corresponds  to  the  average  temperature  of  Si  and  diamond. 

From  Eqs.  (H)  and  (6),  6Si  and  edia  can  be  expressed  as 

eSi  -  (1avTSeav)/(292*a2)  (,2) 

and 


6dia  =  {Tdeav/4000,!)  t(1/a)  "  (1/R)]  (13) 

The  average  temperature  for  the  Si  chip  and  the  diamond  heat  sink  can  be 
obtained  from  the  following  approximation: 

Tav  *  <Tmax  *  <'«> 

For  Si,  is  located  at  the  p/n  Junction  and  Tmin  is  located  at  the 

metal/Si  interface.  For  diamond,  Tj^  is  located  at  the  metal /diamond 
interface  and  ^min  is  located  at  the  diamond/Cu  interface.  Note  that  to 
calculate  Tav,  one  must  know  the  interface  temperatures. 

An  easy  way  to  obtain  the  interface  temperatures  is  to  use  an 
equivalent  circuit.  This  is  Justified  because  heat  conduction  is  similar 
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to  electrical  conduction.  Note  that  Eq.  (1)  is  the  thermal  analog  of  Ohm's 
law.  Both  isothermal  and  equipotential  lines  satisfy  Laplace's  equation. 
The  thermal  resistance  in  degrees  centigrade  per  watt  is  related  by  a 
constant  to  the  electrical  resistance  in  ohms  for  a  given  geometrical 
configuration.  The  temperature  is  equivalent  to  the  voltage,  while  the 
heat  flow  is  analogous  to  the  current  flow.  Figure  5  shows  the  equivalent 
circuit  of  the  heat  flow  of  an  IKPATT  diode  structure  toward  the  heat  sink, 
i.e.,  downward;  this  figure  is  the  electrical  representation  of  Eq.  (3). 

The  bottom  temperature  of  Fig.  5  is  denoted  by  Ta,  the  ambient  temperature 
(T  =  293  K),  while  the  top  temperature  is  denoted  by  T„  «•  AT  (this  is  also 

a 

the  Junction  temperature  Tj  of  the  IMPATT  diode  and  Tmax  of  the  Si  chip) . 
Thus 

Tj  =  Ta  +  AT  (15) 

The  interface  temperatures  can  be  obtained  by  using  the  equivalent 
circuit  of  Fig.  5  as  a  voltage-  or  temperature-dividing  network.  The 
values  of  eSi  and  edia  [Eqs.  (12)  and  (13),  respectively]  can  be  obtained 
by  using  an  iterative  procedure.  From  Fig.  5  note  that  6dia  depends 
upon  e<m  and  vice  versa.  Therefore,  both  values  must  be  calculated 
simultaneously. 

When  the  values  of  8ffietal,  6Cu,  eSi,  and  edia  from  Eqs.  (5),  (7), 

(12),  and  (13),  respectively,  are  entered  into  Eq.  (3),  67,  the  thermal 
resistance  of  the  IMPATT  structure,  becomes 

9T  *  (1a»TSeav)/(292*a2)  *  I,  Vki 

n=  1 

♦  (Tjea,/1*000*)  [( 1/a)  -  (1/R]  ♦  1/(4FkCu)  (16) 

Equation  (16)  represents  the  proposed  thermal  model  of  the  IMPATT 
diode. 
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Fig.  5.  Equivalent  Circuit  of  the  Heat  Flow  of  an  IMPATT  Diode  Structure 
Toward  the  Heat  Sink  (Downward) 


V.  APPLICATION  OF  THE  THERMAL  MODEL 


The  most  important  application  of  the  proposed  thermal  model  of  the 
double-drift  IMPATT  diode  is  in  evaluating  accelerated  life-test  data.  The 
objective  of  accelerated  life  testing  is  to  determine  the  reliability  of 
IMPATT  diodes  by  accelerating  thermally  activated  failure  modes.  Life-test 
data  are  used  to  estimate  the  expected  life  of  an  IMPATT  diode  under  normal 
operating  conditions.  The  thermal  model  as  expressed  by  Eq.  (16)  is  used 
to  define  Junction  temperature,  Si /metal  interface  temperatures,  and  so  on, 
for  three  double-drift  IMPATT  diode  lots:  DMX  109,  DMX  114  and  DMX  117. 
Important  conclusions  are  drawn  from  the  calculations. 

Table  2  lists  the  diode  parameters  of  lots  DMX  109,  DMX  114,  and  DMX 
117;  these  parameters  are  used  to  calculate  the  accelerated  life-test 
conditions. 


Table  2.  Diode  Parameters  of  Lots  DMX  109,  DMX  114,  and  DMX  117 


Diode  Lot  No. 

DMX  109 

DMX  114 

DMX  117 

Junction  Depth,  pm  (Xpn) 

1.17 

0.7 

0.72 

Doping  Concentration,  cm"^ 

na 

5.6  x  1016 

9.5  *  1016 

1.2  x  1017 

nd 

7.4  *  1016 

9.7  *  1016 

9  *  1016 

Chip  Thickness,  pm 

15.5 

18.2 

13.4 

1chip»  vm 

14.33 

17.5 

12.68 

Chip  Radius,  pm 

28.2 

25.4 

25.5 

Breakdown  Voltage  at  20°C,  V 

28.59 

22.63 

22.60 

Space-Charge  Resistance,  Q 

6.51 

4.06 

4.18 

The  accelerated  life  test  is  performed  at  elevated  temperatures  under 
bias  conditions,  where  the  bias  is  used  to  generate  the  heat  necessary  to 
obtain  the  targeted  AT.  Standard  life-test  conditions  are  determined  by 
the  following  method:  (1)  determine  the  thermal  resistance  at  room 
temperature,  0q,  for  all  diode  samples  (room  temperature  is  assumed  to  be 
20°C);  (2)  reverse-bias  the  sample  diodes  so  that 


AT0  =  eoIV 


(17) 


is  the  same  for  all  diodes;  (3)  start  with  ATq  =  305°C,  then  after  a 
predetermined  time  (-168  hr)  successively  increase  the  aTq  to  330,  355,  and 
380°C  or  until  failure  occurs;  and  (4)  record  failure  times  and  temper¬ 
atures  for  all  diodes. 

The  life-test  condition  expressed  by  Eq.  (17)  uses  an  artificial 
Junction  temperature  that  is  denoted  by  TQ.  T 0  can  be  expressed  as 

T.  =  T  ♦  8nl V  =  20  ♦  enIV  (18) 

0  a  o  u 

T q  values  used  for  life  tests  are  325,  350,  375,  and  400°C.  Actual 
Junction  temperatures  are  substantially  higher  than  those  expressed  by  Tq. 

Failure  in  a  life  test  usually  occurs  at  the  metal/Si  interface  either 
at  the  top  or  bottom  of  the  IMPATT  die.  It  is  important  to  know  the  actual 
interface  temperature  where  failure  occurs.  Note  that  because  of  the 
temperature  drop  across  the  Si  layer  between  the  Junction  and  the  inter¬ 
face,  the  temperature  at  the  metal/Si  interface  is  significantly  below  that 
at  the  IMPATT  diode  Junction. 

The  thermal  model  as  expressed  by  Eq.  (16),  together  with  the 
equivalent  circuit  (Fig.  5),  will  be  used  to  calculate  the  interface 
temperatures  toward  the  heat  sink  of  an  IMPATT  structure  under  life  test. 


To  calculate  the  metal/Si  interface  temperature  at  the  top  of  the 
IMPATT  die,  the  equivalent  circuit  of  the  upward  heat-conduction  path  must 
be  used.  The  thermal  resistance  of  the  upward  heat-conduction  path,  eUp, 
can  be  expressed  as 


e  =  e  e  ^  ♦  e  ...  +  e  .  (19) 

up  chip  metal  ribbon  ring 

where  echip  represents  the  upward  thermal  resistance  of  the  Si 
diode,  6metai  is  the  thermal  resistance  of  the  metal  layer,  6ribbon  is  the 
thermal  resistance  of  the  gold  bonding  ribbon,  and  8ring  is  the  thermal 
resistance  of  the  ceramic  or  quartz  ring  used  in  the  assembly.  When 
calculating  the  value  of  8chip,  one  should  use  Eq.  (M)  but  should  replace 
lav  with  lcb|p-  emetal  can  be  ca^ cula ted  after  the  data  presented  in  Table 
1  are  entered  into  Eq.  (5). 

The  equivalent  circuit  of  Eq.  (19)  is  shown  in  Fig.  6;  note  the 
difference  in  the  temperature  values  in  comparison  with  those  in  Fig.  5. 

The  bottom  temperature  is  denoted  by  Ta,  +  AT,.  This  is  also  the  Junction 
temperature  of  the  IMPATT  diode,  which  is  the  same  as  that  expressed  in  Eq. 
(15).  The  top  temperature  is  denoted  by  Ta,.  We  begin  by  examining  the 
expression  for  T j : 


Tj  =  Ta1  ♦  AT,  (20) 

Comparing  Eq.  (20)  to  (15),  we  can  see  that  the  ambient  temperature  and  the 
temperature  rise  (i.e.,  Ta,  and  AT,,  respectively)  are  different.  The 
assembly  of  IMPATT  diodes  explains  this  difference. 

Figure  7  shows  a  detailed  cross  section  of  the  assembly  of  an  IMPATT 
structure.  In  assembly  a  ceramic  or  quartz  ring  is  bonded  on  top  of  the  Cu 
base.  A  gold  ribbon  that  provides  the  top  contact  of  the  diode  is  then 
bonded  to  the  top  of  the  package  ring.  The  inside  radius  of  the  ring  used 
in  the  assembly  is  330  urn,  while  the  outside  radius  is  *<32  urn.  The  radius 
of  the  diamond  heat  sink  is  380  um.  Therefore,  in  assembly  the  ring 
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.  Equivalent  Circuit  of  the  Heat  Flow  of  an  IMPATT  Diode  Structure 
Through  the  Chip,  the  Connecting  Ribbon,  and  the  Ring  (Upward) 
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overlaps  the  diamond  heat  sink  and  the  Cu  base,  as  shown  in  Fig.  7.  As  a 
result,  the  bottom  temperature  at  the  inside  of  the  ring  is  higher  than  the 
bottom  temperature  at  the  outside  of  the  ring.  The  ambient  temperature  of 
the  upward  path,  Ta1,  is  the  temperature  at  the  bottom  of  the  ceramic 
ring.  This  temperature  is  difficult  to  define,  to  say  the  least,  and  can 
only  be  approximated.  The  approximation  used  defines  the  ambient  tem¬ 
perature  Tai  as  the  temperature  at  the  interface  of  the  diamond  heat  sink 
and  the  Cu  base,  instead  of  room  temperature.  Because  of  this  difference  in 
the  ambient  temperature,  the  temperature  rise  must  be  adjusted 
accordingly,  since  there  is  no  change  in  the  actual  Junction  temperature. 

Table  3  lists  the  pertinent  dimensions  of  the  three  diode  lots  and  the 
diamond  heat  sink  used  to  calculate  the  thermal  resistance  values. 


Table  3.  Diode  Junction  and  Diamond  Heat  Sink  Dimensions 
(All  Measurements  are  in  ym  ) 


Lot  No. 

Chip 

Radius 

1av 

■*chip 

Diamond 

Radius 

Heat  Sink 

Height 

DMX  109 

28.2 

1.17 

14.33 

380 

250 

DMX  114 

25.4 

0.7 

17.5 

380 

250 

DMX  117 

25.5 

0.72 

12.68 

380 

250 

Table  4  gives  the  component  values  and  the  total  thermal  resistance  at 
room  temperature  (eQ)  for  lots  DMX  109,  DMX  114,  and  DMX  117.  The  values 
of  and  edia  were  calculated  by  entering  data  from  Table  3  into  Eqs. 

(12)  and  (13).  The  value  of  8ms;tal  was  calculate  with  Eq.  (5)  and  the 
data  shown  in  Table  1.  The  value  of  6qu  was  calculated  from  kCu  =  3.77 
W/cm  K. 
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Table  Component  Values  and  Total  Thermal  Resistance  of  Lots 
DMX  109,  DMX  114,  and  DMX  117 


®Si 

0 

metal 

6dia 

®Cu 

eo 

DMX  109 

0.01604  TSeav 

4.75 

0.02612  Tdeay 

1.75 

14.41 

DMX  114 

0.01183  TSeav 

5.84 

0.02924  Tdeav 

1.75 

15.34 

DMX  117 

0.02071  TSeav 

5.84 

0.02911  Tdeav 

1.75 

15.37 

The  total  thermal  resistance  of  an  IMPATT  structure  at  elevated 
temperatures  is  obtained  by  calculating  the  appropriate  values  of  8^ 
and  6dia.  An  iterative  procedure  is  used  to  define  the  corresponding 
values  of  TSeav  and  Tdeav  when  eSi  and  6dia  are  calculated.  These  values, 
along  with  the  other  thermal  resistance  values  stated  in  Table  4,  define 
the  total  thermal  resistance  eT  of  the  IMPATT  diode  structure;  this 
resistance  corresponds  to  a  junction  temperature.  Note  that  when  an  IMPATT 
diode  operates  or  is  life  tested,  then  AT,  the  rise  in  the  Junction 
temperature,  is  given  by 

AT  =  eTiv  (21) 

where  I  and  V  are  the  diode's  bias  current  and  voltage,  respectively.  The 
value  of  8q  as  given  in  Table  4  is  used  to  define  the  accelerated  life-test 
conditions  expressed  by  Eq.  (17).  Note  that  if  one  replaces  with  e0  in 
Eq.  (21),  then  ATq  is  obtained  instead  of  AT.  This  defines  the 
relationship  between  the  actual  Junction  temperature  and  the  life-test 
conditions. 

To  calculate  the  upward  thermal  resistance  of  an  IMPATT  diode 
structure,  the  value  of  O^bon  and  ering  must  be  obtained.  In  our  case 
the  values  given  by  the  manufacturer  were  used:  6ribbon  =  305.5  or 
299.2°C/W  for  the  ribbon  and  8ring  =  21.8eC/W  for  the  ceramic  package  ring. 


6gi  was  calculated  by  entering  the  value  of  lchip  from  Table  3  into  Eq. 
(12),  and  8metal  was  calculated  from  data  in  Table  1.  Table  5  gives  the 
component  values  and  the  total  upward  thermal  resistance  at  room 
temperature  (e^)  of  lots  DMX  109,  DMX  114,  and  DMX  117. 


Table  5.  Component  Values  and  Total  Upward  Thermal  Resistance  of 
Lots  DMX  109,  DMX  114,  and  DMX  117  (in  ®C/W) 


8  u- 

chip 

0 

metal 

0ribbon 

e  . 
ring 

60up 

DMX  109 

TSeav 

1.75 

305.5 

21.8 

364.2 

DMX  114 

°-2’6  TSeav 

1.75 

299.2 

21.8 

375.7 

DMX  117 

°-213  TSeav 

1.75 

305.5 

21.8 

367.1 

To  obtain  the  total  upward  thermal  resistance  of  the  IMPATT  structure 
at  life-test  temperatures,  the  formerly  mentioned  iterative  procedure  is 
used  to  calculate  TSeav  in  order  to  obtain  e^.  Note  that  when  doing  the 
iterative  procedure  one  must  use  the  ambient  temperature  Tal;  thus  the 
interface  temperature  of  the  diamond  heat  sink  and  the  copper  base  must  be 
determined  beforehand.  This  value  of  8gj  and  the  other  thermal  resistance 
values  given  in  Table  5  are  then  used  to  obtain  the  total  thermal 
resistance  eup. 

By  comparing  the  value  of  eQ  of  Table  4  with  the  value  of  e^p  of 
Table  5,  we  can  see  that  eup  is  indeed  much  larger  than  edown;  thus  the 
upward  heat-conducting  path  accounts  for  less  than  5 %  of  the  total  heat 
conduction.  This  Justifies  the  simplifying  assumption  made  before. 

After  the  thermal  resistance  values  were  obtained,  the  temperatures  at 
five  different  points  of  the  IMPATT  diode  structure  were  calculated  by 
means  of  the  electrical  circuit  analog  (Figs.  5  and  6).  The  five 
temperatures  were  at  (1)  the  Junction  (Tj ) ;  (2)  the  bottom  of  the  Si  chip, 
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at  the  Si  and  metal  interface  (Tbot);  (3)  the  top  of  the  diamond  heat  sink, 
at  the  diamond  and  metal  interface  (Tdia);  (*0  the  interface  of  the  diamond 
heat  sink  and  the  copper  base  (TCu  s  Tfi1);  and  (5)  the  top  of  the  Si  chip, 
at  the  Si  and  the  metal  interface  (Tfc0p).  The  calculations  were  performed 
from  Tq  =  100°C  to  Tq  =  400°C  for  each  lot.  Note  that  the  low  temperatures 
of  the  calculated  values  correspond  to  the  normal  operating  temperatures  of 
IMPATT  diodes,  while  the  high  temperatures  correspond  to  the  life- test 
conditions. 

Figure  8  shows  a  plot  of  the  calculated  Tj,  T^,  Tdia,  TCu,  and  Ttop 
as  a  function  of  Tq  for  Lot  DMX  109.  Both  scales,  the  vertical  and  the 
horizontal,  are  in  °C.  The  plot  assumes  an  ambient  temperature  of  20°C. 
Figure  8  represents  the  thermal  model  of  DMX  109.  Figures  9  and  10 
represent  similar  plots  for  Lots  DMX  114  and  DMX  117,  respectively. 
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VI.  DISCUSSION 


From  Figs.  8,  9.  and  10  it  can  be  seen  that  under  life-test  conditions 
the  Junction  temperature  of  IMPATT  diodes  ranges  from  about  500°C  at  TQ  = 
325°C  to  about  700°C  at  Tq  =  M00°C.  From  the  quoted  junction  temperatures 
it  can  easily  be  concluded  that  IMPATT  diodes  are  subject  to  severe  thermal 
stress  during  life  test. 

In  IMPATT  diodes  life-test  failures  usually  occur  at  the  metal/Si 
interface,  where  temperatures  are  significantly  below  the  Junction  tem¬ 
perature.  Failure  analysis  performed  on  diodes  that  failed  life  test 
indicates  that  the  dominant  failure  mode  is  a  metal/Si  encroachment  through 
the  Junction,  usually  near  the  center  of  the  device.  In  life  test  the 
metal/Si  interface  temperatures  range  from  a  low  of  about  370°C  to  a  high 
of  greater  than  500°C.  These  temperatures  are  high,  exceeding  the  Si-Au 
eutectic  temperature  (377°C)  used  as  a  benchmark.  The  metal/Si  interface 
temperatures  dictate  the  use  of  device  metal ization  capable  of  withstanding 
temperatures  greater  than  370°C. 

The  metalization  of  an  IMPATT  diode  consists  of  a  chromium/platinum/ 
gold  (Cr/Pt/Au)  layer,  where  the  Cr  layer  is  in  contact  with  the  Si 
surface.  There  are  three  binary  interfaces  (Cr/Si,  Pt/Cr,  and  Au/Pt)  whose 
properties  at  elevated  temperatures  affect  the  operation  of  the  device. 

The  Cr/Si  phase  diagram  applicable  to  the  Cr/Si  interface  indicates 
that  Cr  and  Si  may  form  several  binary  compounds,  among  them  Cr^Si,  CrSi, 
and  CrSig.^’®  There  is  little  formation  of  these  compounds  at  temperatures 
below  600‘C.  In  addition,  there  is  little  evidence  that  Cr  enters  into 
solid  solution  with  Si.  This  means  that  Cr  forms  a  good  barrier  in  the 
temperature  range  of  interest  for  IMPATT  diodes  and  is  not  the  cause  of 
thermal  failures. 

Cr  is  also  in  contact  with  Pt,  and  from  the  Cr/Pt  phase  diagram  it  can 
be  seen  that  Pt  dissolves  Cr  through  the  formation  of  CrPt  and  CrPtj.  Of 
these  two  compounds,  the  latter  is  the  more  stable.  Usually  CrPt  converts 
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to  CrPtj  at  temperatures  greater  than  400°C.  Also,  CrPt  and  CrPt^  are 
dissolved  by  Pt.  Table  1  shows  the  usual  thicknesses  of  the  various  metal 
layers.  It  can  be  seen  from  this  table  that,  during  lifetest,  Pt  dissolves 
the  Cr  layer  and  eventually  the  original  Cr/Si  interface  converts  to  a 
Pt/Si  interface. 

The  Pt/Si  phase  diagram  indicates  that  Pt  starts  forming  silicides 
with  Si  at  low  temperatures  (-200°C),  but  accelerated  silicide  formation 
occurs  at  temperatures  above  400°C.^  Furthermore,  Pt  diffuses  rapidly  into 
Si  when  temperatures  exceed  450  to  460°C.^  From  Figs.  8,  9,  and  10,  note 
that  temperatures  exceeding  450°C  exist  in  the  bulk  of  the  Si  even  at  the 
lowest  stress  temperatures,  i.e.,  where  Tq  =  325 °C.  Therefore,  after  Cr 
goes  into  solid  solution  with  Pt,  Pt  starts  to  diffuse  into  Si.  Note  also 
that  since  the  Junction  of  the  diode  is  the  hottest  part  of  the  IMPATT 
structure,  Pt  will  congregate  to  the  p/n  Junction  and  cause  device 
failures. 

The  solubility  of  Au  into  Pt  or  vice  versa  is  negligible  for 
temperatures  under  600°C.  Thus  Pt  forms  a  good  barrier  against  Au.  The 
Importance  of  this  is  that  while  Si  and  Au  form  a  eutectic  at  377°C,  the  Pt 
layer  forestalls  the  Si/Au  eutectic  formation  process  even  though  the 
metal/Si  interface  temperatures  exceed  the  Si/Au  eutectic  temperature. 
Frequently  Pt  is  replaced  by  palladium  < Pd ) ,  but  the  Cr/Pd  and  Si/Pd  phase 
diagrams  indicate  a  great  similarity  between  Pd  and  Pt. 

As  said  before,  the  objective  of  life  testing  IMPATT  diodes  is  to 
determine  the  long-term  reliability  of  the  devices  by  accelerating  ther¬ 
mally  activated  failure  modes.  The  failure  mode  described  above  is  a 
thermally  activated  one.  However,  if  valid  conclusions  are  to  be  reached, 
three  conditions  must  be  met:  (1)  The  failure  mode  that  dominates  during 
life  testing  must  be  the  dominant  failure  mode  during  normal  operation. 

(2)  The  same  failure  mode  must  be  the  dominant  one  for  every  lot.  (3)  For 
valid  projections  to  be  made,  the  life  expectancy  must  vary  in  a  known  way. 
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In  the  case  of  thermal  stress  one  follows  the  Arrhenius  rule,  which 
says  that  the  expected  mean  time  to  failure  (MTTF),  t,  can  be  expressed  as 


t  =  tq  exp(Eft/kT)  (22) 

where  tq  is  a  constant,  Eft  is  the  activation  energy  of  the  thermally 
induced  failure  mechanism,  k  is  Boltzman's  constant,  and  T  is  the  absolute 
temperature.  The  values  of  and  tq  are  determined  by  accelerated  life 
testing.  This  satisfies  condition  (3),  above. 

One  can  see  from  Figs.  8,  9,  and  10  that  during  operation  or  during 
life  test  the  temperatures  at  the  top  and  bottom  metal /Si  interface  are  not 
the  same.  In  addition,  there  is  an  inconsistency  about  the  location  of  the 
hotter  metal/Si  interface;  in  the  case  of  lots  DMX  109  and  DMX  117  it  is 
the  top  contact,  while  in  the  case  of  lot  DMX  114  it  is  the  bottom 
contact.  Since  failure  originates  at  the  location  of  the  hotter  metal/Si 
interface,  thermal  failures  may  iginate  at  either  the  top  or  bottom 
contact.  This  violates  condition  (2),  since  it  can  be  argued  that  lots  DMX 
109  and  DMX  114  are  dominated  by  different  failure  modes.  For  condition 
(2)  to  be  satisfied,  the  hotter  interface  should  always  be  at  the  same 
location,  say,  at  the  bottom  contact. 

Using  the  equivalent  circuits  shown  in  Figs.  5  and  6,  one  can  easily 
calculate  the  top  and  bottom  interface  temperatures.  A  quick  analysis  of 
the  problem  indicates  that  if  the  following  condition  is  satisfied, 


echip/60up  <  0Si/6O 


(23) 


then  the  hotter  interface  will  always  be  located  at  the  bottom  of  the 
die.  In  Eq.  (23)  6Chip  and  represent  room-temperature  values.  We  can 
reduce  Eq.  (23)  to  a  simple  rule  of  thumb  by  requiring  that 


^chip  >  20  ^av 


(24) 
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Equation  (24)  advocates  the  use  of  a  slightly  thicker  diode  chip  than 
is  preferred  today,  if  negative  consequences  can  be  minimized.  The  reason 
for  thinning  the  wafer  is  to  reduce  the  series  resistance  of  the  n+ 
substrate,  thereby  improving  the  conversion  efficiency  of  the  diode. 

Adding  extra  thickness  to  the  wafer,  however,  may  not  reduce  the  conversion 
efficiency  drastically,  because  the  requirement  expressed  by  Eq.  (24)  would 
onlj  increase  the  series  resistance  of  the  IMPATT  diode  from  about  5%  to 
less  than  8J  of  the  space-charge  resistance. 

The  actual  Junction  temperature  Tj  of  a  device  during  life  test 
exceeds  Tq,  the  temperature  defined  by  the  life-test  conditions.  The 
reason  for  this  is  that  the  room-temperature  values  of  and  are 

used  to  define  Tq.  However,  6^,1  anc*  e<jia  are  a  *'unction  of  temperature  and 
increase  with  increasing  temperature.  Therefore,  the  thermal  resistance  of 
the  IMPATT  structure  at  the  life-test  temperatures,  0T,  is  larger  than  6q, 
resulting  in  Tj  >  TQ.  To  reduce  the  difference  between  Tj  and  Tq,  the 
values  of  eSi  and  edia  must  be  kept  small. 

In  addition  to  being  a  function  of  temperature,  0^  is  also  a  function 
of  the  physical  dimensions  of  the  diode  chip.  It  is  directly  proportional 
to  lav,  the  thickness  of  the  Si  layer  between  the  Junction  and  the  heat 
sink,  and  it  is  inversely  proportional  to  a,  the  diameter  of  the  mesa 
diode.  To  obtain  as  small  a  value  of  0^  as  possible,  the  diode  must  have 
a  minimum  l_v,  and  a  maximum  a. 

d  v 

When  operating  as  a  microwave  source,  the  IMPATT  diode  exhibits  a 
negative  ac  resistance;  i.e.,  there  is  a  phase  shift  of  n  between  the 
voltage  and  current  waveforms.  A  phase  shift  of  x/2  is  created  by  the 
avalanching  Junction  and  a  phase  shift  of  x/2  is  created  by  the  drift  of 
the  carriers.  The  width  of  the  drift  region  is  an  inverse  function  cf  the 
frequency;  i.e.,  it  becomes  smaller  with  increasing  frequency.  In 
double-drift  IMPATT  diodes,  the  drift  region  extends  to  both  sides  of  the 
p/n  Junction;  as  a  result  there  is  a  relatively  thick  layer  of  Si  between 
the  Junction  and  the  heat  sink. 
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Since  1&V  is  determined  by  the  thickness  of  the  epitaxial  layer,  it  is 
constant  for  every  diode  fabricated  from  a  given  Si  wafer.  While  1_„  is 
the  same  for  every  diode  in  a  lot,  the  radius  a  of  individual  devices 
within  a  lot  varies  considerably  as  a  result  of  a  final  trim-etch  pro¬ 
cessing  step.  This  means  that  for  a  diode  lot  the  radius  a  of  a  diode  chip 
becomes  an  important  parameter  when  calculating  e^.  Figure  11  shows  a 
normalized  plot  of  0si/TSeav  (normalized  to  lay  =  1  ym)  as  a  function  of 
a.  Note  that  the  plot  of  Fig.  11  is  independent  of  temperature.  Since  the 
radius  of  a  diode  chip  is  usually  20  to  30  urn,  a  small  variation  in  a 
results  in  a  large  variation  in  e^. 

It  can  be  shown  that  the  maximum  radius  of  an  IMPATT  diode  chip,  amax , 
is  a  function  of  the  heat  sink  used  and  can  be  expressed  as 


amax  =  ^  ATmaxkhs^VopJmin^  ^5) 

where  ATmax  is  the  maximum  rise  of  Tj  allowed,  khs  is  the  thermal 
conductivity  of  the  heat  sink,  VQp  is  the  operating  voltage,  and  Jmin  is 
the  minimum  current  density  needed  for  acceptable  conversion  efficiency. 
Equation  (25)  demonstrates  the  need  for  a  heat  sink  to  have  excellent 
thermal  conductivity  and  a  large  AT^jj  (large  Tj). 

The  use  of  a  diamond/copper  heat  sink  has  the  effect  of  increasing  the 
thermal  conductivity  of  the  heat  sink.  To  calculate  the  thermal  resistance 
of  a  heat  sink,  an  expression  for  thermal  spreading  resistance  i s  used. 

The  thermal  spreading  resistance  of  a  heat  sink  is  a  function  of  the  radius 
of  the  diode  chip  and  the  thermal  conductivity  of  the  heat  sink.  Figure  12 
shows  the  difference  between  the  thermal  spreading  resistance  of  a  copper 
or  a  diamond /copper  heat  sink.  The  solid  line  shows  the  thermal  spreading 
resistance  of  a  copper  heat  sink  as  a  function  of  the  radius  of  the  IMPATT 
die,  while  the  dotted  lines  show  the  combined  thermal  spreading  resistance 
of  a  diamond /copper  heat  sink  as  a  function  of  the  die  radius  at  a  diamond 
temperature  of  373  K  (100°C)  and  573  K  (300°C),  respectively,  using  a 
diamond  chip  having  a  radius  of  380  urn.  The  373  K  diamond  temperature 
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corresponds  to  standard  operating  temperature,  while  the  573  K  diamond 
temperature  occurs  during  life  test.  Figure  12  demonstrates  the  advantage 
of  using  diamond  heat  sinks. 

Finally,  the  implementation  of  a  minimum  lav  and  a  maximum  a  leads  to 
the  lowest  possible  junction  temperature  and  the  lowest  metal/Si  interface 
temperatures . 


VII.  CONCLUSIONS 


He  have  developed  a  thermal  model  of  double-drift  IMPATT  diodes  on  a 
diamond  heat  sink  by  using  empirical  formulas  [Eqs.  (10)  and  (11)]  that 
closely  approximate  the  thermal  conductivity  of  Si  and  diamond  (Type  II). 

An  application  of  this  thermal  model  to  three  IMPATT  diode  lots  under 
life  test  indicates  that  extreme  Junction  temperatures ,  in  excess  of  700°C, 
are  reached;  these  temperatures  are  higher  than  those  predicted  by  former 
thermal  models.  In  the  three  diode  lots  modeled,  the  temperatures  at  the 
metal/Si  interface  exceeded  the  benchmark  Au/Si  .eutectic  temperature 
(377°C).  It  was  found  that  at  these  elevated  temperatures  the  diffusion  of 
Pt  into  Si  is  the  cause  of  thermal  failures. 

The  existence  of  such  extreme  Junction  temperatures  during  life  test 
is  the  result  of  the  method  by  which  the  life-test  conditions  are  de¬ 
fined.  Under  life  test  the  IMPATT  diode  is  subject  to  extreme  stress.  In 
order  to  reduce  the  temperature  stress  to  a  more  acceptable  value,  the 
IMPATT  die  must  have  a  minimum  lflV  and  a  maximum  a.  Unfortunately,  there 
are  limits  on  both  of  these  parameters. 

Our  thermal  model  also  indicates  a  need  to  reevaluate  the  life-test 
results  of  IMPATT  diodes.  Prior  thermal  models  that  were  used  to  determine 
T  when  calculating  t0  and  EA  underestimated  the  actual  Junction  temperature 
during  life  test.  Therefore,  by  using  the  thermal  model  developed  here, 
one  should  correct  the  life  expectancy  of  IMPATT  diodes  to  reflect  these 
changes. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  "archi tect -eng ineer ”  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific,  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination ,  thermal  and  structural 
control;  high  temperature  thermomechanics,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  excimer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -f ield-of  -view  rejection, 
applied  laser  spectroscopy,  laser  chemistry,  laser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verif ication,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 

Electronics  Research  Laboratory:  Microelectronics ,  solid-state  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communications; 
microwave  semiconductor  devices,  raicrowave/mill imeter  wave  measurements, 
diagnostics  and  radiometry,  microwave/milliraeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments. 

Space  Sciences  Laboratory:  Magnetospher ic,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  Infrared  astronomy, 
infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  Che  earth’s  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
i ns t rumen tat  ion. 


